Measuring radial flow of partonic and hadronic phases in relativistic heavy ion 

collision. 
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It has been shown that the thermal photon and the lepton pair spectra can be used to estimate 
the radial velocity of different phases of the matter formed in nuclear collisions at ultra-relativistic 
energies. We observe a non-monotonic variation of the flow velocity with invariant mass of the 
lepton pair which is indicative of two different thermal dilepton sources at early and late stage of the 
dynamically evolving system. We also show that the study of radial velocity through electromagnetic 
probes may shed light on the nature of the phase transition from hadrons to QGP. 

PACS numbers: 25.75.-q,25.75.Dw,24.85.+p 



I. INTRODUCTION 

The numerical simulation of QCD equation of 
state (EoS) predicts that nuclear matter at high 
density and/or temperature are composed of quarks 
and gluons due to asymptotic freedom and screen- 
ing of colour charges [L, 2, 3]. Enormous experimen- 
tal efforts have been made to produce such a par- 
tonic state of a matter, called Quark Gluon Plasma 
(QGP) by colliding nuclei at ultra-relativistic ener- 
gies. Careful theoretical investigations have been 
performed to understand the existing experimental 
data [4|and predictions for the forthcoming experi- 
ments [5| have also been made. 

The hot and dense matter formed in the partonic 
phase after ultra-relativistic heavy ion collisions ex- 
pands in space and time due to high internal pres- 
sure. Consequently the system cools and reverts 
back to hadronic matter from the partonic phase. 
Initially (when the thermal system is just born) the 
entire energy of the system is thermal in nature and 
with progress of time some part of the thermal en- 
ergy gets converted to the collective (flow) energy. 
In other words during the expansion stage the total 
energy of the system is shared by the thermal as well 
as collective degrees of freedom. The evolution of the 
collectivity within the system is sensitive to the EoS. 
Therefore, the study of the collectivity in the system 
formed after nuclear collisions will be useful to shed 
light on the EoS [1, 0, [1] of the strongly interacting 
system at high temperatures and densities. 

It is well known that the average magnitude of 
radial flow can be extracted from the transverse 
momentum (pt) spectra of the hadrons. However, 
hadrons being strongly interacting objects can bring 
the information of the state of the system when it 
is too dilute to support collectivity. On the other 
hand electromagnetic (EM) probes i.e. photons and 
dileptons are produced and emitted from each space 
time points. Therefore, estimating radial flow from 



the EM probes will shed light on the time evolution 
of the collectivity in the system. This is demon- 
strated by NA60 collaboration through dilep- 
ton measurements in In-(-In collisions at SPS energy. 
The slope of the transverse mass spectrum of lepton 
pairs, T e ff of invariant mass M can be related to 
the space-timed averaged quantities like radial flow 
velocity u r and temperature T av as T e g ~ T av + Mvf. 
T c f j estimated from dilepton spectra M shows a dif- 
ferent kind of behaviour 0, E EE E d as 
compared to that from hadronic spectra. The ef- 
fective temperature extracted from transverse mass 
spectra of dileptons increases linearly with invariant 
mass M up to p-peak and then falls (the PHENIX 
data does not show this trend [II). In a recent 
work [17] we have shown that the ratio (i? cm ) of 
the px spectra of photons to lepton pairs has an 
advantage over the individual spectra because some 
of the uncertainties or model dependence pertaining 
to the individual spectra gets canceled in the ratio. 
Hence the ratio can be used as an efficient tool to 
understand the state of an expanding system. In the 
present work we focus on the extraction of the radial 
flow from .Rem- We also argue that the simultaneous 
measurements of photons and dileptons will enable 
us to estimate the value of v r for various invariant 
mass windows of the lepton pairs. The v r obtained 
from the analysis of both the spectra vary with M 
non-monotonically. Such a behaviour may be inter- 
preted as due to the presence of two different kinds 
of thermal sources of lepton pairs of the expanding 
system. 

The paper is organized as follows. In section II 
the ratio of thermal photon and dilepton productions 
has been discussed. In section III the evolution dy- 
namics of the hot fireball system with specific initial 
conditions and EoS is outlined. The discussions in 
sections II and III will be very brief as the details are 
available elsewhere |17[ • The results are presented in 
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section IV. Finally section V is devoted to summary 
and discussions. 



II. ELECTROMAGNETIC PROBES 

The ratio, R em of the px spectra of thermal pho- 
tons to dileptons can be written as follows [13] : 
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The numerator (denominator) is the invariant mo- 
mentum distribution of the thermal photons (lepton 
pairs). In Eq. [Upt, V and M denote the transverse 
momentum, rapidity and the invariant mass of the 
lepton pair. The summation in Eq. [1] runs over all 
phases through which the system passes during the 
expansion. (d 2 R/d 2 pTdy)i and {d 2 R/d 2 p T dydM 2 )i 
are the static rates of photon and diiepton pro- 
ductions from the phase i, which is convoluted 
over the expansion dynamics through the space- 
time integration over d A x. The integration over 
M is done by selecting invariant mass windows - 
M m i n < M < M max appropriately and we define 
< M >= (M mm + M max )/2. 

The rate of thermal diiepton production per unit 
space-time volume per unit four momentum volume 
is given by [H M 0, E3 



dR 
~¥p~ 
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(2) 



where a is the EM coupling constant, ItoII^ is the 
imaginary part of the retarded photon self energy 
and fBE{E,T) is the thermal phase space factor for 



Bosons. L(p 2 )=(l + ^S-Vl 



4^V arises from the 
p 2 

final state leptonic current involving Dirac spinors of 
mass m. The real photon production rate can be ob- 
tained from the diiepton emission rate by replacing 
the product of EM vertex 7* — > , the term in- 
volving final state leptonic current and the square of 
the (virtual) photon propagator by the polarization 
sum for the real photon. For an expanding system E 
should be replaced by u^p^ where p^ and are the 
four momentum and the four velocity respectively. 



A. Thermal photons 

The photon production rate has been evaluate d by 
various authors [12] using Hard Thermal Loops [23| 
approximations. The complete calculation of emis- 
sion rate of photons from QGP to order 0(a, a s ) has 
been done by resuming ladder diagrams in the effec- 
tive theory [2J]. This rate of production has been 



considered in the present work. A set of hadronic 
reactions with all isospin combinations have been 
considered for the production of photons [25|, [2(| H3| 
from hadronic matter. The effect of hadronic dipole 
form factors has been taken into account in the 
present work. We have checked that the high pt 
(~ 2 — 3 GeV) part of the thermal photon spectra 
is dominated by the contributions from QGP phase 
for large initial temperature. 



B. Thermal dileptons 

The lowest order process producing lepton pairs 
is q and q annihilation. For a finite temperature 
QCD plasma, the correction of order a s a 2 to the 
lowest order rate of diiepton production has been 
calculated in [28|, [2t|, which is considered in the 
present work. For the low M diiepton production 
from the hadronic phase we consider the decay of 
light vector mesons p, cj and (A as considered in [l7| . 
The continuum part of the vector mesons spectral 
functions have been included in the present work [30L 

It is well known that the contributions from the 
QGP phase dominates the M spectra of the lepton 
pairs below p-peak and above the </>-peak if no ther- 
mal effects of the spectral functions of the vector 
mesons (see [3(1 HH, [33| for review) are considered. 



III. EVOLUTION DYNAMICS 

In the collision of two energetic heavy ions a large 
amount of energy is dumped into a small volume. 
The space-time evolution of the matter has been 
studied using ideal relativistic hydrodynamics (34| 
with longitudinal boost invariance [35j and cylindri- 
cal symmetry. The initial energy density (e(ri,r)) 
and radial velocity (v(ri, r)) profiles are same as our 
earlier studies [l7| . The value of transition temper- 
ature (T c ) is taken as 192 MeV as obtained in lat- 
tice QCD calculations [36| . Alth oug h a much lower 
value of T c is also predicted in 37]. However, we 
have found that the dependence of R em on T c is 
weak. In a first order phase transition scenario we 
use Bag EoS for the QGP phase and for the hadronic 
phase all the resonances with mass < 2.5 GeV have 
been considered [Hj]. 

In the present work we have considered the ini- 
tial and freeze out conditions that reproduced the 
hadrons [3{|, the photon [i(| and diiepton spectra 
for RHIC energy [4l|. The values of initial ther- 
malization time, 7j = 0.2 fm/c, initial tempera- 
ture Ti — 400 MeV and the freeze-out temperature 
Tf = 120 MeV have been taken as the input to the 
calculation. For LHC we have taken Ti — 700 MeV, 



3 



> 



Z 

X! 



10" 



10" 



10" 



10 



10" 




1 1.5 2 
P T (GeV) 



FIG. 1: The Pt spectra of photons and dileptons from 
hadronic and quark matter at RHIC energy. The dilep- 
ton spectra is obtained by doing M integration from 
M =1.0 GeV to 1.4 GeV 



Ti = 0.08 fm/c which gives the hadron multiplicity 
dN/dy = 2100 @. 



IV. RESULTS AND DISCUSSION 

In Fig. Q] the photon and dilepton spectra have 
been displayed for RHIC conditions. Results indi- 
cate that the photon spectra from QGP dominates 
over its hadronic counterpart for pr > 1.5 GeV. The 
dilepton from QGP and hadrons are comparable in 
magnitude for entire range of pr for M ~ 1.2 GeV 
(this is because of the inclusion of the continuum of 
the vector meson spectral functions [30j, [31| without 
the continuum the quark matter part dominates). 
However, for M ~ 0.75 GeV the dileptons from the 
hadronic matter are overwhelmingly large compared 
to quark matter contributions (not shown in the fig- 
ure). Therefore, an appropriate selection of px and 
M will be very useful to characterize a particular 
phase of the system. 

Now we consider the variation of the ratio, R em 
as a function of px for different invariant mass win- 
dows. The results are shown in the Fig. [2] and Fig. [3] 
below. The variation of R cm with respect to px can 
be parametrized as follows: 

R em = A 3 [^} B3 exp[C 3 {m T - PT )} 

PT 



(3) 

where A3, B3 and C3 are constants and Mt, the 
transverse mass of the lepton pair is defined as, 
Mt — y/p T + M 2 . It is observed that the ra- 
tio decreases sharply and reaches a plateau beyond 



Pt > 1.5 GeV. This behaviour of R cm as a func- 
tion of pt can be understood as follows: (i) for 
Pt » M Mt ~ Pt and consequently R em ~ A3 
giving rise to a plateau at large pt- The height 
of the plateau is sensitive to the initial tempera- 
ture of the system [17|. (ii) For pr < M R cm ~ 
exp{—pT /T c f f) /p T 3 indicating a decrease of the ra- 
tio with pt (at low pt) as observed in the Fig. [5] and 

Fig. m 

To indicate the effect of the radial flow velocity, v T 
we have evaluated the R cm with and without radial 
flow (see Fig. [5] and Fig. [3]). In case of vanishing 
radial flow the ratio can be parametrized as follows: 



R\ m = A 1 [^} B ^exp[C 1 (m T - PT ) 

PT 



(4) 



Here C\ contains the information of the average tem- 
perature, T av of the system. 

In case of vanishing radial flow velocity the inverse 
slope of the photon and dilepton spectra represent 
the average temperature, T av of the system. How- 
ever, in case of non-zero radial flow the inverse slope 
contains the effect of average temperature as well as 
that of v r . Therefore, the difference in the slopes of 
the two cases will enable us to estimate the amount 
of collectivity in the system. 

As mentioned before for large initial temperature 
transverse momentum distribution of photons from 
QGP dominates over its hadronic counterpart for 
Pt > 1.5 GeV. However, in case of dileptons one 
has to select both the M and the pt windows to 
observe QGP. For example the thermal dileptons 
from hadrons dominate over those from QGP for 
M ~ 0.75 GeV. Therefore, for estimating the radial 
velocity in the hadronic phase we chose pt ~ 0.5 
GeV and M ~ 0.75 GeV for demonstrative pur- 
pose. Similarly a pt and M windows may be se- 
lected where contributions from QGP dominates. 

The exponential slope of the ratio (C3) can be 
related to the individual slopes of photons CT e //l ) 
and dileptons {T~ f \ 2 ) as follows: 



C 3 x (m T - pt) 



m,T 



Pt 



T e ff2 T e ff! 

writing the effective (blue shifted) temperatures of 
the photon spectra and dilepton spectra as 



T e ffl - T av \ 



'(l + O 



. (1-Vr)' 
Teff2 = T av + Mv r 2 



(5) 



we obtain, 
C 3 x(m T - Pt) 



m,T 



Pt 



T av +Mv r 2 T av ^/(l + v r )/(l~v r ) 
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FIG. 2: R cm as a function of pr with and without ra- 
dial flow for invariant mass 0.6 < M(GeV) < 0.9. The 
spectra with radial flow is normalized to the one without 
radial flow at pr — 0.5 GeV 
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FIG. 3: -Rom as a function of pr as in the above figure 
for invariant mass 1.0 < M(GeV) < 1.4. 



Further simplification leads to 



bT av + c = 



(6) 



where a, b and c are functions of v T . Solving Eq.[S]for 
a given C3 , M and px we obtain v T as a function of 
the average temperature. The results are displayed 
in Figs. |31 and [5] for initial conditions of RHIC and 
LHC energies for invariant mass and px windows 
indicated. The contributions in the M and px win- 
dows shown in Fig. [4] are dominated by the hadronic 
phase i.e. from temperature range T c ~ 192 MeV to 
Tp ~ 120 MeV. The radial velocity increases sharply 
with decrease in T &v in the hadronic phase. 



FIG. 4: Variation of v t with T av for M = 0.75 GeV 
and pt = 0.5 GeV. The solid (dashed) line indicate the 
results for RHIC (LHC) for EoS with first order phase 
transition. The line with asterisk (dotted line) stands 
for RHIC (LHC)for an EoS which excludes the mixed 
phase. 



We have evaluated v r with a (continuous) EoS 
where the mixed phase does not appear. In this case 
the v r is larger than the one obtained for a strong 
first order phase transition ( Fig. [|]). Indicating that 
the presence of the mixed phase (of hadrons and 
QGP) characterized by zero sound velocity slowed 
down the expansion of the system, resulting in a 
lower radial flow. Therefore, extraction of v T from 
experimental data will be useful to understand the 
nature of the transition. 

In Fig. [5] the radial velocity is displayed for (aver- 
age) temperature range which is dominated by QGP 
phase. The results indicate a moderate v r for RHIC 
but a large v r is achieved even in the QGP phase 
for LHC energies, in fact a fast increase in v r is ob- 
served for T av close to the transition temperature in 
case of LHC. The v r for LHC is much larger than 
RHIC because of the longer life time and larger in- 
ternal pressure of the partonic phase in LHC than 
RHIC. 

In a first order phase transition scenario the QGP 
formed in heavy ion collisions return back to hadrons 
through a first order phase transition. The temper- 
ature changes continuously from Tj to Tp. We esti- 
mate the average values of the radial velocity (vi so th) 
on the constant temperature surfaces determined by 
the conditions: T(r, r) = Tg, for various values Tg. 
The variation of (vi so th) with T5 is depicted in Fig.[S] 
both for RHIC and LHC energies. v- lso th for LHC is 
larger than RHIC because of higher initial tempera- 
ture and hence internal pressure. In contrast to the 
results shown in Figs. [Hand \5\ the variation of Vi SO th 
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FIG. 5: Variation of w r with T a „ for M = 1.2 GeV and 
p T = 0.5 GeV at RHIC and LHC energies for EoS with 
a first order phase transition. 
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FIG. 6: Variation of average radial velocity of the fluid 
on the constant temperature surface. 



with Ts is not measurable as it does not depend on 
the kinematic variables, pr and M. The expansion 
is slower in the hadronic phase because of the softer 
EoS as compared to the QGP phase. For given T c 
and Tp the life time of the hadronic phase is larger 
for softer EoS - allowing the system to develop large 
radial flow as evident from the results depicted in 
Fig. [5] for the low temperature part. The effective 
temperature extracted from the ratio is displayed in 
Fig. [7] as a function of M for RHIC energy. T c f / in- 
creases with M up to the p-peak and then decreases 
beyond p mass. The reduction of T e // beyond p is 
indicating the dominance of the radiation from the 
high temperature phase in the high M region. For 
LHC, however, no clear reduction of T e ff beyond p 



- peak is observed (Fig. [H]). At LHC the average 
temperature and the flow velocity in the early phase 
(from where large M pairs originate) are large (see 
Fig. 2]). Hence the combination of both large v r and 
large T av does not allow T c ff to fall above the p- 
peak. The dependence of individual spectra on Tp 
is quite strong, however, we have observed that the 
slope of the ratio is insensitive to Tp and also to T c . 
The slope of the ratio does not change when the pa- 
rameters like Tp changes from 0.120 GeV to 0.150 
GeV and T c from 0.192 GeV to 0.175 GeV. 

Eliminating T &v from Eq. [5] and taking the val- 
ues of T ffi and T e //2 from photon and dilepton 
spectra one can obtain the variation of func- 
tion of M. The results are shown in Fig. [5] for 
RHIC and LHC energies. A non-monotonic be- 
haviour of v r with M is observed. A similar non- 
monotonic behaviour is observed in the elliptic flow 
(V2) of photons as a function of transverse momen- 
tum [42|, E3|. Comparison of dilepton production 
from QGP and hadronic sources [17| indicate that 
in the low M (< m p ) and high M(> m^) region the 
emission rate from QGP dominates over its hadronic 
counter part if the medium effects on the vector me- 
son spectral functions are neglected. In other words, 
for a dynamically evolving system the low and high 
M pairs are emitted from early QGP phase, whereas 
lepton pairs with M around p - mass are emitted 
from the late hadronic phase. Therefore, low and 
high M domains represent early time where v r is 
low and the M ~ m p domain represent late time - 
where v r is large - giving rise to the observed varia- 
tion in Fig. [9]- indicative of a two different kinds of 
source in early and late times of the evolving sys- 
tem. For M ~ 1.2 GeV the flow velocity is not 
very small since this window is populated by both 
hadronic and partonic contributions almost equally. 
Again at LHC energy the partonic phase life time 
is more which favors the development of larger flow 
compared to RHIC energy. It is important to note at 
this point that for LHC, although the slope C3 does 
not show a clear non-monotonic behaviour with M , 
v r does so. Because as described, before the slope 
C3 depends not only on v r but also on T av and both 
are large in the partonic phase at LHC. 

The two time scales - the life time of the par- 
tonic phase (tqgp) and the time an inward moving 
rarefaction wave takes to hit the centre of the cylin- 
drical geometry decide whether radial flow play any 
important role in the partonic phase or not. The 
later time scale is defined as t to ~ R/c s where R 
is the transverse size of the system and c s is the 
velocity of sound. If tqgp ~ r IW then v x will be 
large in the partonic phase. Therefore, an increase 
in Ti (tqgp (x Ti) will increase the radial flow in the 
partonic phase if the initial and the critical tempera- 
tures are kept fixed. However, an increase in Ti from 
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LHC energies. 
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n to T2 produces same flow if the Ti decreases by 
a factor (T2/T1) 1 / 3 . For a fixed Ti an increase in Ti 
will increase the effective slope as evident from the 
right panel of Fig. |TD] Therefore, the slope of the 
ratio may be used effectively to estimate the value 
of initial thermalization time. 



V. SUMMARY AND DISCUSSIONS 

It has been shown that the pr distribution of ther- 
mal photons and lepton pair spectra may be used 
simultaneously to estimate the magnitude of the ra- 



dial velocity of different phases of the matter formed 
in nuclear collisions at ultra-relativistic energies. Ju- 
dicious choices of the kinematic variables e.g. the 
invariant mass and the transverse momentum win- 
dows may be selected to estimate the flow velocity 
in the partonic and hadronic phases of the evolving 
matter. It has been observed that for RHIC and 
LHC energies the flow velocity increases with invari- 
ant mass up to the p peak beyond which it decreases. 
The T c f f may not decrease with mass beyond p peak 
if the average temperature and the flow velocity are 
large in the partonic phase as in case of LHC energy. 
By doing a simple analysis of photon and dilepton 



7 



spectra we have extracted the radial flow velocity 
for various invariant mass windows. v r varies with 
M non-monotonically. We argue that such a varia- 
tion indicates the presence of two different types of 
thermal sources of lepton pairs. 
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